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Is there a standard measuring rod in the Universe? 
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ABSTRACT 

The CaltechJodrcll Bank very long baseline interferometry (VLBI) Surveys give de- 
tailed 5 GHz VLBI images of several hundred milliarcsecond (mas) radio sources, and 
the full width at half-maximum angular sizes of the corresponding compact cores. 
Using the latter, I have constructed an angular-diameter/redshift diagram comprising 
271 objects, which shows clearly the expected features of such a diagram, without 
redshift binning. Cosmological parameters are derived which are compatible with ex- 
isting consensus values, particularly when the VLBI data are combined with recent 
Baryon Acoustic Oscillations observations; the figures are presented as indications of 
what might be expected of larger samples of similar data. The importance of beaming 
and relativistic motion towards the observer is stressed; a model of the latter indi- 
cates that the emitting material is close to the observers line of sight and moving with 
a velocity which brings it close to the observers rest frame. With respect to linear 
size, these objects compare reasonably well in variance with the absolute luminosity 
of type la supernovae; the efficacy of the latter is improved by the brighter-slower and 
brighter-bluer correlations, and by the inverse-square law. 
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1 INTRODUCTION 

, Cosmological parameters are now known to a remark- 
■ able degree of precision, particularly those derived from 
Wilki nson Microwave Anisotropy Probe fWMAP ) observa- 
! tions (iDunklev et al.1 l2008l : iKomatsu etakl 120081 ) . in com- 
' bina tion with observations of Type la supernovae (SNe 
la) jRiess et al.' 2004; ' Astier et alTl2006l : iRiess et al.l l2007l : 
IWood-Vase~ et al,, .200^ and the imprint of Baryon Ac- 
coust ic Oscillations (BAG) on t he distribution of galax- 
ies (iBlake fc Glazebroo^ |2003| ; ISeo fc EisensteinI l2003l : 
lEisenstein et al.l 20051 ); current values are fim = 0.279 ± 
0.013 and JIa ~ 0.721 ± 0.015 for the matter and vacuum 
parameters respectively (unless otherwise noted, all confi- 
dence limits quoted here are 68 per cent ones) . Whereas the 
WMAP and BAG approaches are of relatively recent origin, 
the magnitude/redshift approach has a long and somewhat 
varied history, for want of a class of objects with similar 
absolue magnitudes (see e.g. [Weinberg (|l972l ) for a review 
of early work); however, over the last decade the latter ap- 
proach has achieved spectacular success, with the discovery 
of SNe la as accurate standard candles l|Riess et al] 1 19981 ; 
ISchmidt et allll998l ; [Perlmutter et al.lll999l i 

In contrast, the angular-size/redshift approach has had 
a surprisingly modest impact, despite significant efi'orts in 



this direction, using extra-galactic radio sources in their 
several guises as putative standard measuring rods. Early 
work co nsidered classical double radio sources as suitable 
objects (|Legg| Il970l ; iMilevI ll97lU Kellermann' 'l972') , which 
approach continues to this day (Daly 1994; Buchalto r et al.l 
ll998l ; [Dalv et al.l2007^ . Here I will re-examine ultra-compact 
radio sources as standard measuring rods, with angular di- 
ameters in the milliarcsecond (mas) range, and linear sizes 
of order several parsecs. Their advantages in this context 
were first highlighted by iKellermannI (|l993h ; these objects 
are much smaller than their parent active galactic nuclei, so 
that their local environments should be similar and reason- 
ably stable, at least over an appropriate redshift range. They 
are energetic and shortlived, with central engines which are 
reasonably standard objects (black holes with masses close 
to 1.5 X IO'^'^Mq). In these respects they have much in com- 
mon with SNe la, albeit w i th life times of centuries rather 
than months. IKellermannI l|l993l ) presented angular sizes 
for a sample of 79 milliarcsecond (mas) sources, obtained 
using very long baseline interferometry (VLBI) at 5GHz. 
VLBI imag es tyically show a compact core surrounded by 
debris, and IKellermannI (|l993h defined a characteristic an- 
gular size as the distance between the core and the most 
distant component having a peak brightness greater than 
or equal to 2 per cent of that of the core. Typical lin- 
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ear sizes are pcQ and Kellermann noted that the 

corresponding angular-size/redshift diagram is compatible 
with the once-favoured flat cold dark matter (CDM) model, 
f2m = 1, Ha = 0. Ifowever, for a critique see Pearson et 
al. (1994); these authors present a similar sample (with 
very little overlap) and use the same measure of angular 
size, which shows virtually no change over the same redshift 
ranee. iGurvits. Kellermann fc FrevI (|l99^) have examined a 
much larger 5 GHz sample (330 objects), with results which 
do not resolve this conflict. I suspect that the discord here 
is due to the definition of angular size, which is too sensitive 
to t he detail s of so urce structure. 

iGurvitsI (l994) presented a large VLB I compilation, 
based upon a 2.29 GHz survey undertaken bv lPreston et al.l 
l|l985l ): the compilation lists a rough measure of angular size 
based upon fringe visibility; this measure should be less sen- 
sitive to the details of source structure, and more repre- 
sentative of the compact core, which is usually the domi- 
nant component wit h respect to radio luminosity (see e.g. 
iHenstock et~aLlll995h . Gurvits considered a subset compris- 
ing 258 sources divided into 12 redshift bins over the range 
0.501 ^ 2 ^ 3.787, and found marginal support for a low- 
density CDM cosmological model, assuming that selection 
and evolutionary effects can be ignored. Gurvits considered 
only models with SI a = 0. Using ex actly the same data 
set (k indly supplied by Dr. Gurvits), Ijackson fc DodgsonI 
|l993) extended the analysis to the full f2m- f^A plane; the 
situation is very degenerate with respect to choice of Q.a, 
with both 95 per cent and 68 per cent confidence regions ex- 
tending well into the region SIa < 0. Nevertheless, marginal- 
izing over Oa > gives 0.11 ^ Sim ^5 0.54, providing clear 
support for a low-density CDM model. Imposing spatial flat- 
ness allowed a much more definitive statement to be made: 
0-1 <^m< 0.3, later refined t o = 0-24 + 0.09/ - 0.07 
(95 per cent confidence limits) l|jacksonll2004l ). Fl 



2 DATA AND RESULTS 

The object here is to consider a more precise definition of 
angular size. The Caltech-Jodrell Bank fiat-spectrum (CJF) 
sample is a complete 5 GHz fiux-density-limited sample of 
293 flat-spe ctrum sources, com plete according to the follow- 
ing criteria (iTavlor et al.ll 19961 ): 

(i) Flux density at 4850 MHz 54850 > 350 mjy. 

(ii) Spectral index a ^ —0.5 (S oc frequency"). 

(iii) Declination 5 ^ 35° (1950 coordinates). 



I Ijackson h Jannettal ll200d) . 



(iv) Galactic latitude ^ 10° 



1 Ho = lOO/i km sec-i Mpc"! 
^ Tighter figures are presented 
where a larger sample i s exa mined (613 objects), produced by 
updating iPreston et al.l lll985l) with respect to redshift, and find- 
ing some of the missing flux densities from elsewhere; this work 
places significant constraints on Qrn and Qa without assuming 
flatness. However, in retrospect we have doubts about the suit- 
abilty of the extra data, and about th e efficacy of re dshift binning, 
see discussion later; the precision in Ijackson fc Jannetta (200a ) 
has been overstated. 



Figure 1. Angular-diameter /redshift diagram for 271 sources 
from the composite PR+CJ1-I-CJ2 sample. The cyan line cor- 
responds to rim = 0.236, Ha = 1 - Hm = 0.764 and d = 2.94/i-l 
pc, see text. 



iTavlor et all (|l996D present a VLBI image for each 
successfully imaged source, and give accurate full width 
at half-maximum angular major and minor axes (a and 
6) for each component therein, with the compact core 
clearly identified; the major axis of the latter will be 
taken as the measure of angular size in this investigation. 
In fact I have used a somewhat larger sample than the 
CJF one. The latt er is, in part, a subset of three earlier 
samples, the PR dPearson fc ReadheadI Il98ll. Il988l) CJl 



(jPolatidis et~allll995l: iThakkar et al.lll995l: 



and CJ2 jTavlor et al.lll994l : 



Henstock et al 



Xu et al 



1995l i samples 



A composite comprising the latter three is nominally com- 
plete with respect to crite ria (i), (iii) and (iv ) above, but 
not the spectral index one. ITavlor et all (|l996l ) selected the 
CJF sample from this composite by imposing the spectral 
index limit, and also added 18 further sources which had 
been missed in the earlier surveys, which additional sources 
meet all of the above criteria. In the interest of a modest 
increase in numbers I have used the f ull compo s ite sa mple 
PR-FCJ1-HCJ2 as list ed in tables 3 o f IXu et al.l (Il995l) and 
iHenstock et al. (Il995l) and tab les 4 of iPearson fc ReadheadI 
(| 19881 ) and ITavlor et al.l (|l994l), plus 17 of the th e 18 further 
sources, as listed in table 3 of Tavlor et al.l (|l996t l , giving 322 
objects in all. Of these 6 are assigned major axes formally 
equal to 0.00 milliarcsec, presumably because they are be- 
low the resolution limit of the VLBI system; the said sources 
have been discarded. I have updated the redshift list using 
the NASA/IPAC Extragalactic Database (NED), and for 
the remaining 316 objects find 271 redshifts. Fig. [l] is a plot 
of major axis a against redshift z for these sources. Despite 
the spread the expected qualitative features are reasonably 
clear, namely a diminishing size from 2 = to z ~ 1, folowed 
by a gradual increase, with a minim um somewhe re between 
z ~ 1 and z 1.5, first predicted bv lHovld (|l959l ). The plot 
in Fig. [1] does not appear to be resolution limited. 

Before turning to quantitative matters, a second dia- 
gram is instructive. Fig. [2] shows linear size plotted against 
radio luminosity for the 271 sources in Fig. [1] taking Sim = 
0.27, Qa ~ 0.73. The cyan points correspond to z < 0.5, the 
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radio luminosity (h W/Hz) 

Figure 2. Plot of linear size against radio luminosity for 271 
sources from the composite PR+CJ1+CJ2 sample; the cyan 
points have z < 0.5, the blue ones z ^ 0.5. 



blue ones to z ^ 0.5. In a flux-limited sample sources ob- 
served at large redshifts are intrinsically the most powerful, 
so that a correlation between linear size and radio lumi- 
nosity would introduce a selection effect. The high-redshift 
population shows no obvious evidence of such a correlation, 
and gives every indication of being statistically stable with 
respect to linear size. In the low-redshift case this is clearly 
not so. Fig. [2] gives a clear illustration of why ultra-compact 
radio sources with z < 0.5 are of no value in this context. 



first n oted by Guryitd (Il994l '). see also Ijackson &: DodgsonI 
([l993) and lJacksonI (|2004 ). 

In what follows I will work with data points which cor- 
respond to individual sources, rather than putting the latter 
into redshift bins. Although binning is popular and can re- 
veal trends in data which are not otherwise obvious to the 
eye, with regard to quantitative statistical analysis there 
is, in principle, no advantage in such a procedure; what is 
gained by having smaller error bars is lost by having fewer 
points. Appearances can be deceptive, particularly if num- 
bers are small, and results can be overly sensitive to the 
choice of bins. I will concentrate on spatially flat ACDM 
models, characterized by Qrn and a characteristic linear size 
d associated with the source population. I will use objects 
in the range 0.5 55 z 55 3.5; the upper limit on z removes 3 
points, one of which is an outlier which has an inordinate 
effect upon the statistical analysis. 

Giving each point equal weight and taking log a as ordi- 
nate, the best-fitting fiat model is f2m = 0.931, flA ~ 0.069 
and d — 2.33h~^ pc, close to the erstwhile canonical model 
Jim = 1, flA ~ 0. However, there is another parameter which 
I believe is an important source discriminator, namely the 
axial ratio r = b/a. This belief is based upon an astro- 
physical model, discussed at length in Jackson (2004). Ac- 
cording to this model the underlying source population con- 
sists of compact symmetric objects ((Wilkinson ct al. 1994), 
comprising central low-luminosity cores straddled by two 
mini-lobes. The compact components which are the basis 
of this study are identified as cases in which the lobes are 
moving relativistically, and are close to the line of sight, 



when Doppler boosting allows just that material which is 
moving towards the observer to be seen. As z increases 
a larger Doppler factor is required; it turns out that the 
latter approximately cancels the cosmological redshift, so 
that the observed component is seen in its rest frame. This 
is a very important effect, because the measured size of 
mas source components is known to increase linearly with 
wavelength (jMarscher fc Shaffeij |l980: Pearson fc ReadheadI 
1981). Without the said effect mas angular-size/redshift dia- 
grams would show something like the so-called Euclidean be- 
haviour, angular size proportional to 1/z (because the emit- 
ted frequency is (1 + z) x the received one), w hich bedev- 
illed early work on classical double radio sour ces (|Legg||l970l : 
lMilevlll97ll : iKellermannll 19721 : |jacksonlll973h . It is quite re- 
markable that Fig.[T]shows no trace of such behaviour, which 
fact is a striking confirmation of the proposed model. 

The ideal image will correspond to a head-on approach, 
and an axial ratio close to unity. I have examined samples 
which include only those objects with r J5 rc, and find that 
the best-fitting value of fim depends upon the cut-off ratio 
rc in a systematic fashion. The said value falls from 0.931 
at Tc = to 0.253 at Vc = 0.35, and thereafter remains 
reasonably constant until we run out of sensible numbers; 
fim = 0.28 ± 0.04 for 0.35 ^ < 0.6, and Q.^ = 0.29 ± 0.08 
for 0.35 ^ 55 0.7. As a representative example I present 
results for the case Vc = 0.4, 0.5 55 z ^ 3.5, comprising 
128 objects, for which the best figures are flm = 0.236 and 
d = 2.94/;,"^ pc; the corresponding curve is shown in Fig. [1] 
A fixed standard deviation a is attached to each point, being 
defined by — residual sum-of-squares/(n — p), where n = 
128 is the number of points and p = 2 is the number of fitted 
parameters; the appropriate value in this case is 

a = 0.252, (1) 

which value is used to calculate values at points in pa- 
rameter space. The mid-grey lines in Fig. |3] show confi- 
dence regions in the Qm- d plane. Marginalizing over d gives 
f7m = 0.24 + 0.40/ -0.15. 

This somewhat indeterminate result is in large part oc- 
casioned by the absence of standard objects in the redshift 
range < z < 0.5, which would otherwise determine a 
model-independent normalization of the characteristic lin- 
ear size d. The deficiency can be remedied by combining 
recent BAO observations with the above data. I have in 
mind measures of th e BAO scales at z = 0.2 and z = 0.35 
(|Percival et all |2007| ). which constrain values of the hy- 
brid distance Dv{z) = [(1 + zfD\{z)cz/H{z)]^^^ , giving 
I>v(0.35)/Dv(0.2) = 1.812 ± 0.060, which ratio does not 
refer directly to the size of the acoustic horizon at recom- 
bination. The joint mas-f-BAO confidence regions are repre- 
sented by the filled blue/cyan areas in Fig. [3] Marginalizing 
over d gives = 0.10 + 0.09/ - 0.06. The light-grey verti- 
cal line is the right-most 95 per cent BAO confidence limit; 
the rest are formally at negative values of f2m, including the 
best-fitting line f2m ~ —0.13. 



3 CODA AND CONCLUSIONS 

Ultra-compact radio sources do comprise a reasonably stan- 
dard unit of linear size, and in conjunction with BAO they 
are close to being useful for cosmological investigations. How 
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Figure 3. Joint mas/BAO confidence regions in tlie f2ni— d plane, 
68 per cent (blue) and 95 per cent (cyan); the mid-grey lines refer 
to mas sources alone; the light-grey vertical line is the right-most 
95 per cent BAO confidence limit. Data relate to 128 sources from 
the composite PR-I-CJ1+CJ2 sample, with z ^ 0.5 and axial ratio 
> 0.4. 



well do they compare with SNe la in this respect? Given a 
radio source at redshift z with linear size d and angular size 
6, its luminosity distance is 



d 



(1 + ^) 



(2) 



As apparent and absolute magnitudes are related by 

m = M -f- 5(log]^Q dL — 1) (d_L in parsecs), (3) 

an equivalent mock supernova is defined by 

m-mas = M -f 5[logi(, d + 2 logio(l + z) ~ logn, 9 - 1]. (4) 

A mock magnitude/redshift diagram is produced by as- 
signing approprate fixed values to M and d, as in 
Ijackson fc Jannettal (|2006l ). Thus a variation 56 at fixed z 
is equivalent to a variation 



5mmas = — 551ogi„6'. 



(5) 



The equivalent dispersion is thus Smmas ~ ±1.26 mag, from 
equation ([T|. For SNe la the observed figure is (5msN '~ 
±0.38 ma g, believed to be larg ely intrinsic rather than pho- 
tometric l|Hamuv et al.|[r996al lb[). The present mas perfor- 
mance thus appears to be inferior by a factor of about 
3. In fact SNe la do rather better than this, because 
their effective dispersion is red uced by the brighter-slower 
and brighter-bluer corr elations (|Phillipsll 19931 : iHamuv et al.l 
Il995l: iGuv et al.l |2005| ). typically to &msn - ±0.20 mag 
llAstier et al.l |2006| ~ and 5msN ~ ±0.27 mag (|Riess et al.l 
I2OO4I . the gold subset). The latter figures are believed to 
comprise photometric and residual intrinsic errors in roughly 
equal proportions. It is instructive to look at the basic physi- 
cal parameters in the two cases: the linear size d and absolute 
luminosity L. The variation in L is given by equation ((3)1 as 



5L 

- = -log 



10 X 0.4 X 5msN, 



the equivalent for mas sources being 



^ = log, 10 X (5(logu)( 



' logg 10 X 0.2 X Sniu 



(6) 



(7) 



At first sight equations ^ and ^ present a puzzle; why, 
for equivalent intrinsic variations 5L/L and 5d/d, is (Jmmas 
larger than (5msN by a factor of 2? The answer is that it 
is the inverse-square law, implicit in equations (|3) and ([4]), 
which puts the angular approach at a disadvantage. The 
numerical values are not dissimilar: 5d/d ~ 0.58 and 5L/L ~ 
0.35; the latter is reduced by a factor of about 1.6 by the 
above-mentioned correlations, whereas the former is in effect 
increased by a factor of 2 by the inverse-square law. 

I have looked at spectral index a as a possible size dis- 
cr iminator for mas sou rces, using the integrated values listed 
m iTavlor et all l| 19961 ). but find that the correlation is too 
weak to be useful, over the range of a values encompassed 
here. The central black hole mass Mbh would be an interest- 
ing parameter in this context, which could be determined 
by the correlation between Mb h and velocity dispersion 
of the the host ellipt ical galaxy i Ferrarese fc Merrittll200"ol : 
iGebhardt et al.|[2000l ). or that between Mbh and the Sersic 
luminosity concentration index of the galaxy llGraham et al] 
I2OOII : iDriver et al1l2006l : iGraham fc Driveij|200'i1 ). The ap- 
propriate information is not currently available; its acqui- 
sition would be a worthwhile undertaking. The variance in 
measured size probably has significant contributions from 
both intrinsic and instrumental effects, so that an increase 
in relolving power would improve matters. There is certainly 
scope for increasing the sample size; the Galtech-Jodrell 
Bank surveys cover only about 20 per cent of the full sky, 
so that a five- fold increase would be possible without chang- 
ing the flux limit. The comoving volume encompassed by a 
complete survey of fiat-spectrum sources with fiux-density 
limit 5"; increases roughly as S^^ , over the redshift range 
considered here. 

The mas/BAO combination is a natural one, in that it 
allows cosmological parameters to be determined by data 
which are local {z < 4) and exclusively angular. Any dis- 
cernible differences between parameters so determined and 
those determined by supernovae would be a manifestation 
of differential selection or evolutionary effects, or of ef- 
fects relating to the astrophysics or possibly the fundamen- 
tal physics of li ght propagati on over cosm ological distances 
l|Bassett fc Kun z 2004a b; B urragell20o3 ). On more general 
grounds a new approach is always of some value, even if its 
weight is relatively low, because its systematic errors be- 
come random ones when the new technique is added to an 
ensemble of existing ones. 

Interested parties can obtain copies of the data set 
used in this investigation by sending an email request to 
j ohn . j ackson® unn .ac.uk. 
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